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Abstract 

The solutions to the Einstein-Klein-Gordon equations without a 
cosmological constant are investigated for an exponential potential in 
a Bianchi VIq metric. There exists a two-parameter family of solutions 
which have a power-law inflationary behaviour when the exponent of 
the potential, k, satisfies k'^ < 2. In addition, there exists a two- 
parameter family of singular solutions for all k‘^ values. A simple 
anisotropic exact solution is found to be stable when 2 < k^. 
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power-law inflation. 
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Inflationary theories claim that the isotropy of the Universe can be ex¬ 
plained by assnming an inflationary expansion in the early nniverse. This 
belief is based on the “Cosmic no-hair” theorems []^. Heusler [Q, however, 
proved that by considering the entire evolntion of the scalar held conpled 
to the gravitational held and assuming a very large class of reasonable po¬ 
tentials, the evolution of the anisotropic Bianchi models was similar to that 
found by Collins and Hawking P]: only those Bianchi models that have as 
particular solutions the FRW models isotropize. 

One of the potentials which has received special attention is that of the 
Liouville form. It appears in the Jordan-Brans-Dicke theory and does not 
belong to the class of potentials analyzed by Heusler. The asymptotic be¬ 
haviour of the Bianchi models with an exponential-potential scalar held has 
been studied recently 0 using the techniques developed in p. When the 
constant k that gives the slope of the potential is less than a/2 then the space- 
time inhates and isotropizes. However, when > 2, the models that can 
possibly isotropize are those of Bianchi types I, V, VH or IX. The isotropiza- 
tion of several scalar held Bianchi models has been recently studied in 0 
and 0. 

In 0 the Bianchi I models with exponential potential were analyzed re¬ 
ducing the problem to one third order diherential equation. In @ a non-local 
transformation was used to linearize this equation and the general solution 
was found. It allows a better insight into the problem and shows a damped 
oscillatory behaviour which leads to an ehective negative cosmological con¬ 
stant. Several singular solutions representing universe models that have ei¬ 
ther anisotropic or isotropic Friedmann-Robertson-Walker hnal stages were 
found. There are also solutions which avoid the initial singularity and oth¬ 
ers with a hnite time span. In this paper we shall extend this analysis to 
the Bianchi VIq type, investigating the behaviour of the solutions near the 
singularity and its asymptotic stability in the far future. 


2 


The metric corresponding to a Bianchi VIq cosmological model can be 
written in the following form: 


+ G{t) {e^dx^ + e ^dy^^ . (1) 

The pertinent Klein-Gordon and Einstein field equations for the metric (|^) 
are as follows: 
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We will consider a homogeneous self-interacting scalar field (0 = 0(t)) with 
an exponential potential V = Ae^'^ . As can be easily shown 
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is a first integral of the Klein-Gordon equation and m is an arbitrary inte¬ 
gration constant. From Eqs.(|-|^), we get 


1 k"^ ^ 

G^G-KGG^-GGG + -GG‘^ + m^G = 0; K= -, (6) 

2 4 2 '^ 

First we will consider the special case in which = 0, and after that, we 
will pay attention to the most general case ^ 0. 


m2 = 0 case 

To analyze the asymptotical behaviour of the solutions we define the new 
variables h = ^ and dr] = hdt = dlnG, so that Eq.(^ now reads 
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h‘^h" + (1 + K)h‘^h' + Kh^ + h) = 0, (7) 

with h' = ^. The simplest solution of equation (|^) is ^ for K > 0. 
Hence, we get 


G = Goev^*, (8) 

while the remaining component of the metric / and the scalar held (j) can be 
obtained from Eqs.(|^) and (|^) respectively. 

Now, we assume that the function G vanishes or diverges at f = 0 and its 
leading term (this assumption is justihed in Ref.|]^), is given by 


G(f)~GoG. (9) 

Substituting this expression into Eq.(|^), we can see that the last two terms 
are negligible compared with the hrst four terms. So 
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is the equation that determines the values of the exponent n. Inserting this 
value of n and Eq.(^ in Eqs.(|^) and (|^) , we obtain the approximate scalar 
held and metric in the limit f —> 0. 

Eor 1 < K which corresponds to values of 6 < with 0 < n < 1, we have 
A < 0 for Eq.(H), so the potential must be negative. However, it is positive 
for K <1. 

For — 7< A' < 0 which corresponds to values of <2, with —oo<n<—2, 
the metric coehcients G and e-^ diverge when f —> 0. 

The next term in the expansion of the function G, can be obtained solving 
Eq. (]^) when the last two terms are neglected. Its general implicit solution is 
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where r = / dt. 

For — I < iF < 0 and 1 < K the function G diverges (vanishes) in the limit 
r —>• dzcx) according to G2 > 0 or C2 < 0. In this case we get 

G{t) ~ at + ( 12 ) 

There exists a two-parameter family of solutions that behaves as G ~ \t\^■ 
Also, G vanishes at a hnite time Tq for iF < 1, provided that G 1 G 2 G 3 > 0. 
Its approximate expansion is 


G ~ at -I- h\t\^ ^. (13) 

Thus a two-parameter family of solutions behaves as G ~ t. 

The behaviour and stability of the solutions can be investigated by writing 
Eq.(|^) as the equation of motion for a dissipative or antidissipative system, 
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Equation (JT^) presents local minima when hg = for iF > 0 (i.e. > 2). 

As the dissipative term given by the right-hand side of Eq.(p^ is negative 
dehnite in the asymptotic regime, the bracket on the l.h.s. of Eq.(^) dehnes 
a Liapunov Function |]^,|^,|]^. Then, the exact solution (^ is stable for 
t ^ 00 and for any initial condition. Studying the behaviour of the solutions 
around these equilibrium points we have: 

For iF > 0 (fc^ > 2) there is a two-parameter family of stable solutions that 
behaves as (§). When iF > 1 it can be shown that the solutions cut the h 
axis in the phase plane {h, h) an inhnite number of times. Therefore they 
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spiral in the phase plane around the constant solution h^. For —\<K<^ 
the solutions do not cut the h axis or they cut it once. 


m? 0 case 


We are going to proceed now in the same manner as we did in the previous 
case. In the limit t 0, Eq.(^ with 7 ^ 0 gives 
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From this equation we see that the term which comes from GG^ in can 
be neglected and the remaining approximate equation can be solved using 


the method described in [15 . In this case, we obtain 
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and 9 = -GJG-^dt. 

/? < 1/4: For —^ < iF < 0 and 1 < iF the coefficient of the metric G diverges 
(vanishes) at a finite time 9 = 9q when sgnCi 7 ^ sgnC' 2 . In this case, we get 


+ for — ^ < iF < 0 and 1 < iF, (19) 

This two-parameter family of solutions behaves as when t 0. Also, 

G vanishes in the limit 6 ^ —> ±cxo for iF < 1 and its expansion is 

2 2 

~ + for iF < iFo = 1 + 7 ^, (20) 
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where = C [1/2 ± X]/{K-1). 

f3 > 1/4: Eq.(p!7D describes damped oscillatory solutions. These solutions 
are compatible only with an effective negative cosmological constant p|. 

In the general case Eq. can be written as the equation of motion for a 
dissipative or antidissipative system as follows 
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where now the “potential” is 


and the local extreme points will be given by 


( 22 ) 


hl = 


1 + 2m^e 
2 ^ ■ 


(23) 


Requiring that K > 0, the extreme points are actually minima and the 
right-hand side of Eq.(|2l|) is negative definite, so that the bracket on the 


l.h.s of Eq. (plD defines a Liapunov function. Near the equilibrium points 
the approximate equation governing the trajectories on the phase plane, for 
f —*• cxo, is given by Eq.(^3D whose solution is 
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showing that the final anisotropic solution is asymptotically stable and co¬ 
incides with the solution (IsD . For > 2 the late-time behaviour of an 


inhomogeneization of a Bianchi I model [[I^ is like that of Bianchi VIq. 

The limit behaviour of the scalar curvature when t 0 has two different 
possibilities depending on the values of K: 



R ~ 



for — i < iL < 0 and 1 < iL, 
for 0 < r: < Kq. 


(26) 
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We conclude that there exists a two-parameter family of singular solutions 
which describe a universe that begins from (or ends in) a singularity for any 
values of the constants m, C and k"^. For these solutions the scalar curvature 
diverges when t > 0. However, for — | < iF < 0 it vanishes asymptotically 
and the solution remains regular. 

We can also see that the shear-expansion ratio, 


a _VQ fG-G 
0 ~ T/G + 2G’ 


(26) 


which is usually considered as an indicator of the isotropization, does not 
vanish for the asymptotical solution given by (|]). Therefore, we can say that 
none of the solutions with k‘^ > 2, for which that study was valid, isotropize. 


Some results of this paper are related to those ones obtained in P, 0 
and |T^. Solutions with k"^ <2 show a power-law inflationary behaviour ||16|| 
while solutions with > 2 do not inflate or isotropize, since Bianchi VI 
models do not have FRW as particular solutions 
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